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Abstract 

Motivated by a model of pseudo-Majoron dark matter, we show how the breaking of a global 
symmetry that acts nontrivially in lepton generation space can lead to a viable pseudo-familon dark 
matter candidate. Unlike the pseudo-Majoron, the pseudo-familon in our model decays primarily 
to charged leptons and can account for the excess observed in the cosmic ray electron and positron 
spectra. 
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I. INTRODUCTION 



The compelling evidence for non-baryonic dark matter has motivated a vast literature on 
possible dark matter candidates. Successful models must lead to the correct relic density, 
as required, for example, by the observed rotation curves of galaxies. However, depending 
on what other phenomena one wishes to explain by the presence of dark matter, the details 
of viable scenarios can differ substantially from model to model. In the present work, 
the phenomena we will attribute to the presence of dark matter is the striking rise in the 
cosmic ray positron fraction, from approximately 7 to 100 GeV observed in the PAMELA 
experiment jl| and the apparent excess in the total electron and positron flux, up to about 
~ 1 TeV, as measured by Fermi-LAT and H.E.S.S. s|. While the discrepancy between 
the observed spectra and predicted backgrounds might have a purely astrophysical origin 
(for example, if nearby pulsars are present j^j, 5|), dark matter annihilation p] or decays 7] 
can also provide a source for the excess electron and positron flux. (A recent review of the 
literature can be found in Ref. In models with dark matter that is not exactly stable, 
fits to the data suggest that the dark matter must decay primarily to leptons and have a 
lifetime of C(10 26 ) s. 

Among the many interesting proposals for dark matter candidates in the recent literature 
is one in which a pseudo-Majoron is identified as the dark matter candidate js]. The Majoron 
is a goldstone boson of spontaneously broken lepton number and the pseudo-Majoron is 
the pseudo-goldstone boson that is obtained when the same spontaneously broken global 
symmetry is also explicitly broken by a small amount. In the model of Ref. the explicit 
breaking is accomplished via soft terms in the scalar potential, and the sense in which this 
breaking is small is that the mass scale of these terms is much smaller than the scale at which 
lepton number is spontaneously broken. Although the Majoron can decay to standard model 
neutrinos, the lifetime is extremely long, due to the high-scale associated with the breaking 
of the global symmetry. While we will argue that the soft symmetry-breaking terms assumed 
in the original pseudo-Majoron proposal do not lead to a viable model, we also show how 
this difficulty can be overcome. Hence, the Majoron can serve as a dark matter candidate 
and perhaps also provide detectable signals in neutrino telescopes. 

The model of Ref. [9|, however, is not one which can account for the cosmic ray electron 
and positron excess. Here, we will show that it is possible to construct a similar model 
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Fields: L\ e R l > 2 e R 3 v R <t>H 



Charges: +4 



4 +4 



-4010 



TABLE I: Charge assignments under the U(l) global symmetry. The superscripts indicated gen- 
eration number. 

in which the dark matter candidate is also the pseudo-goldstone boson of a spontaneously 
broken approximate global symmetry, but one which leads to couplings of the dark matter 
candidate to charged leptons rather than neutrinos. Unlike lepton number (which treats all 
generations of leptons identically) our global symmetry will not. Hence, the dark matter 
candidate that we propose here is technically a pseudo-familon, at least as far as the lepton 
sector is concerned. The model we discuss in the next sections provides an existence proof 
that pseudo-familon dark matter is a viable possibility that can accommodate the cosmic 
ray data. 

II. THE MODEL 

We assume the existence of global symmetry G = U(l) that acts nontrivially in lepton 
generation space. The charge assignments of the fields are summarized in Table HI Notice 
that the global symmetry distinguishes charged leptons of the third generation from those of 
the other two. The complex scalar field has charge +1 and its vacuum expectation value 
(vev) completely breaks the symmetry 



The standard model Higgs doublet H is neutral under G and plays its conventional role in 
electroweak symmetry breaking and the generation of fermion masses. 

The spontaneous breaking of the global symmetry leads to a goldstone boson ip, which 
we shall identify via the conventional nonlinear field redefinition 



that softly break G, the familon (p becomes massive and can annihilate via couplings that 
are not suppressed by the scale v^, regions of the model parameter space exist where the 
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familon is sufficiently long lived and has the necessary relic density to be a viable dark 
matter candidate. 

Before proceeding to evaluate the model that we have just defined, it is worth clarifying 
some of the underlying assumptions. First, we assume the existence of a global symmetry. 
While such symmetries can be violated by generic quantum gravitational effects, it is also 
true that they can arise accidentally in low-energy effective theories as a consequence of the 
gauge symmetries in an ultraviolet completion. Hence, as in the Majoron model of Ref. 9[, 
it is reasonable to consider the consequences of relatively simple global symmetries below the 
Planck scale that lead to interesting phenomenology. Secondly, we employ the term familon 
since our global symmetry G acts nontrivially in lepton generation space. This does not 
imply, however, that G is intended to provide a complete understanding of fermion Yukawa 
textures by itself. Instead, we require that G not restrict the lepton Yukawa textures too 
severely, so that the present construction is compatible with a wide range of possible flavor 
models; these may introduce additional symmetries in generation space. Finally, the model 
we present is non-supersymmetric, as a matter of convenience. Possible solutions to the 
problem of fine-tuning in theories with fundamental scalars are well known, but are largely 
independent of the issues that we consider here. Depending on what we learn from the Large 
Hadron Collider, one can modify the model accordingly. 

The symmetry G partly restricts the lepton Yukawa couplings. The charged lepton 
Yukawa matrix Yl, defined by 



- £ D L L Y L He R + h.c. 
indicates that Yl must have the global charges 
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The entries with non-vanishing charges arise via higher- dimension operators, so that at 
lowest order 
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where the and fry are coefficients that are at most 0(1), and A is the lepton flavor scale, 
the scale at which the higher- dimension operators are generated. This form clearly can 
parametrically accommodate the charged lepton Yukawa couplings. However, one can do 
better: the breaking of additional symmetries can provide an explanation for the charged 
lepton mass hierarchy. For example, consider a Z 6 flavor symmetry under which the three 

1...3, 



generations of right-handed charged leptons transform as ej^ — > u l+2 e^ , for % 
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Assuming (<f) F ) ~ (4>) ~ A 2 A, where A 
viable charged lepton Yukawa texture 
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0.2 is of order the Cabibbo angle, one obtains the 
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and the ip£ + £ coupling 
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where £ = (e, fj,, t) t . Note that the vanishing entries in Eq. ( 12. 8ft are not corrected at higher 
order. Since the diagonal blocks of Eq. ( 12.7ft are neutral under G, any higher-order correc- 
tions to these entries that involve <fi must do so via the product ^</>, which is independent 
of tp. Note also that the matrices in Eqs. ( 12. 7p and ( 12. 8 j) are not proportional to each other 
and cannot be diagonalized simultaneously. Due to the extreme smallness of the 13 and 23 
rotation angles required to diagonalize Y F , one finds that the largest entries of the familon 
coupling matrix are the same as in Eq. ( 12. 81) . though generally with different order-one 
coefficients. In the mass eigenstate basis, one finds 
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where we have only retained order-one coefficients Bij for the non-negligible terms. The (p 
decays are lepton flavor violating; the leading decay channels are tp — > e~r + , e + r~, yU~r + 
and yU + r~. The dark matter lifetime of (9(10 26 ) s is obtained when the operator coefficient 
satisfies 
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which implies A ~ 10 18 GeV, up to order one uncertainty. We will take A ~ M*, the reduced 
Planck mass, henceforth. One should keep in mind, however, that the scale A could arise via 
renormalizable physics, for example, the integrating-out of vector-like fermions with masses 
just below M*. 

In the pseudo-Majoron dark matter model, the Majoron decay amplitude is proportional 
to neutrino masses, which accounts for the longevity of the dark matter candidate. One 
might naively expect that a generic familon in the charged lepton sector would decay too 
quickly, due to the much larger charged lepton masses, to be a viable dark matter candi- 
date. The crucial feature of the model we have presented is that our familon couplings are 
proportional to off-diagonal elements of the charged lepton mass matrix. These elements 
can be arbitrarily small, leading to the long dark matter lifetime required by the cosmic ray 
data. 

Neutrino masses in the model arise via the see-saw mechanism. The Dirac and Majorana 
mass terms are defined by 



- C D v c R M RR u R + (L L Y LR Hu R 
The entries of Ylr have the global charges 
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and hence arise all at the same order (0) 4 /A 4 ~ A 8 ; we may therefore write Y^ R = \ s Yl R , 
where Y^ R is an arbitrary matrix with order-one entries. The matrix M RR is not restricted 
at all by G. Hence, the flavor structure of the neutrino sector is largely unconstrained and 
can be chosen to accommodate the data. On the other hand, the overall light neutrino mass 
scale, m u , follows from the seesaw formula 

\ 16 (H) 2 
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where Mr characterizes the right-handed neutrino mass scale. Choosing m 



0.05 eV 



(which is compatible with the data on atmospheric neutrino oscillations PLOj), one finds 
Mr ~ (9(10) TeV. This scale is high enough so that the see-saw formula is an accurate ap- 
proximation and the heavy, right-handed neutrino mass eigenstates present no phenomeno- 
logical problems. 

III. SCALAR COUPLINGS 

Let us first consider the U(l) invariant portion of the scalar sector of the theory. Imagine 
we employ the linear decomposition <fi = (v^ + (f) r + i \/2, where 0j is the massless degree 
of freedom. Let H represent the standard model Higgs doublet, and h the Higgs boson. In 



the pseudo-Majoron model of Ref. [9|, it was argued that the coupling $<j)H^H leads to 
interactions of the form h 2 <frf and vh<pf, where v/y/2 = (H). Since these are the standard 



Higgs portal couplings for a scalar dark matter candidate, the authors of Ref. js| argued 
that the appropriate relic density could readily be achieved. 

This argument, however, is not correct. Using the nonlinear redefinition of 0, Eq. (12.21) . 
one notices that the U(l) invariant (jy<j) is independent of the goldstone mode tp; hence, 
one concludes using this representation that there are no couplings of (f to h that are 
unsuppressed by the scale v^. Both the linear and non-linear representations of the goldstone 
boson, however, should lead to the same physical results. This puzzle can be resolved in 
the linear theory by carefully evaluating the scalar couplings in the mass eigenstate basis. 
For example, one finds couplings of the form vh<f)f and v^tfir 0|, where <p r is (mostly) a 
super-heavy scalar state with mass of order v^. One cannot simply neglect the coupling to 
<p r : there is mixing between h and <p T such that h ~ h and <p r ~ 0{v/v^) ho, where h 
is the Higgs mass eigenstate. Thus, one finds two contributions to the h <f>? vertex, each 
with a coefficient of 0(v). It is straightforward to verify that they cancel, up to corrections 
of 0(v 3 /vl). The absence of couplings that are unsuppressed by the high scale implies 
that one cannot achieve a sufficient annihilation cross section via the couplings in the U(l) 
invariant sector of the theory. 



The previous observations suggest that the desired Higgs- 
the soft-symmetry breaking sector of the theory. As in Ref. 



mortal couplings must arise in 
9[], we assume that the global 



symmetry is an approximate one due to the existence of renormalizable interactions with 
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mass dimension three or less that explicitly break the symmetry. Such interactions must 
provide for an adequate dark matter mass and annihilation cross section, while not leading 
to rapid dark matter decays. To obtain the dark matter mass scale that is suggested by 
the cosmic-ray data, we take the new dimensionful couplings to be TeV-scale in size. The 
familon ip is nonetheless correctly categorized as a pseudo-goldstone boson since the ratio of 
the scales of explicit to spontaneous symmetry breaking is small, (10 3 GeV)/(10 17 GeV) ~ 
10~ 14 . The simplest way to allow for unsuppressed couplings that lead to dark matter 
annihilation is to introduce a real scalar singlet So and the CP-invariant symmetry-breaking 
terms m\(f) 2 + m 2 Sq <fi 2 + h.c. One may consistently impose a Z 2 symmetry that eliminates 
soft terms that are odd in 0, while allowing all the other Lagrangian terms that we have 
considered thus far. Note that m 1 can always be made real by a global U(l) rotation, 
while m 2 is taken real as a parameter choice; this satisfies 't Hooft's criterion for technical 
naturalness since the scalar sector then has an enhanced symmetry, namely a Z 2 where 
ip — y —ip. It follows that the scalar-sector interactions involve only even powers of pseudo- 
goldstone boson field and therefore do not introduce an avenue for rapid decays . The term 
proportional to m 2 includes the interaction rri2S ip 2 /2, where S represents the fluctuation of 
the singlet about its vev. Provided that the necessary tunings of dimensionful couplings are 
employed so that the scalar S and the Higgs field h are light, then the low-energy effective 
theory will contain the terms 

V D \m%^ + \m 2 s S 2 + l -m 2 S V 2 + ^Sh 2 + • • • , (3.1) 

where the dimensionful parameters shown are TeV-scale in size. Generically, S and h will 
mix, so that the familon can annihilate to standard model particles via the s-channel ex- 
changes of both mass eigenstates. A complete study of the allowed parameters space that 
yields the correct relic density is beyond the scope of his letter (a detailed study of a similar 



sector can be found in Ref. Here we will focus on a well motivated limit: For a Higgs 

boson that is very similar to that of the standard model, the h-S mixing should be small. 
If the singlet is an order of magnitude heavier than the Higgs (i.e., similar in mass to the 
familon), then we would expect this to be a reasonable approximation. With mixing effects 



1 Of course, one may obtain scalar decay channels via diagrams involving the extremely small ip£ + £~ 
couplings which violate this Z-2- However, these decays can be neglected since their branching fractions 
are loop suppressed relative to the leading modes discussed earlier. 



S 



neglected, one sees that the interactions in Eq. (13. ip nonetheless lead to the dark matter 
annihilation channel (ftp — > hh, via s-channel exchanges of the singlet S. The annihilation 
cross section is given by 
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where the singlet width, for ms < 2 m™, is given by 
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For example, for the parameter choices m v = 2 TeV, mj = 3.5 TeV, m2 
// = 2.0 TeV and rrt/j = 125 GeV, we find that the freeze out condition 
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2.5 TeV, 

(3.4) 



gives x = m<p/Tf ~ 28, where n^® is the equilibrium number density, (av) is the thermally 
averaged annihilation cross section and Tf is the freeze out temperature; we find that this 
leads in turn to the present dark matter density Vt£>h 2 ~ 0.1, as desired 2 . This example 
demonstrates that there are regions of the full parameters space of the model where the 
appropriate dark matter relic density can easily be obtained. 



IV. COSMIC RAYS 

The branching fractions for the leading familon decay channels are determined by the 
couplings £>i3 and _B 2 3 in Eq. (I2.9p . i.e., 

C D 4iV2 ^<^A 18 [B 13 eP R r + B 2 ^P R r] + h.c. , (4.1) 
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where Pr = (1 + 7 5 )/2, and are given by 

B(e~T + ) = B(e + T~) = 0.5 
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B(^r + ) = BQj+t-) = 0.5 , - ■ ■ ,„ (4.2) 



Since the total must sum to one, we can parameterize Eq. ( 14. 2ft in terms of the branching 
fraction to ^t+: B(^t + ) = B(fi + r~) = f and B(e~T+) = B(e + r~) = 1/2 - /, with 



2 This analysis is quite standard; a more detailed discussion can be found in Ref. 
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FIG. 1: Positron fraction and total electron plus positron flux for the parameter choice / = 0.375. 
The best fit predicted spectrum (x 2 /d.o.f = 0.61) is shown, corresponding to the familon mass and 
lifetime n%n = 2 TeV and r ip = 2 x 10 26 s, respectively. Data from PAMELA Fermi LAT and 
H.E.S.S. |3j are also shown. The dashed lines indicate the H.E.S.S. band of systematic uncertainty. 



< / < 1/2. Hence the dark matter mass m v , lifetime t v and the value of the parameter / 
determines the contribution from dark matter decays to the cosmic ray spectra. We compute 
the energy distribution of cosmic ray electrons or positrons using PYTHIA 6.4 13J]; when 
more than one decay channel is possible 



dN e ± 
dE 



* 3 



i -J 



dN e ± 
dE 



where (dN e ± / dE),^ is the energy spectrum obtained given the primary decay <p — > 
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with 



the subscript indicating lepton flavor. The energy spectrum in Eq. (14. 3p is then input into a 
standard diffusive model for propagation of electrons and positrons through the interstellar 
medium. This analysis is now standard to the literature on decaying dark matter and the 



detai 
Ref. 



s, including our assumptions on background fluxes, are identical to those found in 



141 ] ; we refer the reader to this reference for a detailed discussion. An example of 



typical results that can be obtained in the model are shown in Fig. [U Here we have fixed 
the parameter / at 0.375 and have performed a least-squares fit to determine the optimal 
values of and t v . Note that the "step" in the predicted total flux around 1 TeV is due 
to the admixture of e ± r =F in the primary decay, which contributes electrons and positrons 
of a fixed injection energy, « m^/2, to the spectrum. In the example shown, this feature 
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is within the total experimental uncertainty of the data in this region, and can be made 
smaller for larger values of /. 

V. CONCLUSIONS 

Pseudo-Majoron dark matter is an interesting proposal that cannot account for the 
anomalies in the cosmic ray electron and positron spectra. We have shown that the breaking 
of a global symmetry that acts nontrivially in lepton generation space can lead to a viable 
dark matter candidate that does decay primarily to charged leptons. The long lifetime of 
the pseudo-familon dark matter candidate in the model we have proposed is related to the 
smallness of off-diagonal entries in the charged lepton Yukawa matrices; no phenomenolog- 
ical consideration prevents these entries from being small, an outcome that is realized due 
to the symmetries, charge assignments and the choice of mass scales in our model. We 
have argued that unsuppressed Higgs-portal couplings can arise via soft symmetry-breaking 
terms involving a gauge singlet scalar, so that the correct pseudo-familon relic density can be 
achieved. We have shown also that the predicted cosmic ray spectra can fit the current 
data from the Fermi-LAT, PAMELA and H.E.S.S experiments. The model we have pre- 
sented serves as a proof of principle; a more elegant implementation within a comprehensive 
framework for the origin of lepton flavor is a direction worthy of investigation. 
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